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Overcasts  seen  from  below  seldom  are  uniform,  unchanging  cloud  shields,  yet  little  is  known  about  their 
short-term  photometric  variability  (periods  <2  h).  Visible-wavelength  spectra  of  daytime  and  twilight 
overcast  skies  measured  at  30-s  intervals  reveal  unexpected  temporal  variability  in  horizontal  illumi¬ 
nance  Ev  and  zenith  luminance  Lv.  Fourier  analysis  of  these  time  series  shows  peak  fluctuations  at 
periods  of  2-40  min.  Factors  such  as  cloud  type  and  optical  depth,  presence  of  fog  or  snow,  and  instrument 
field  of  view  can  affect  overcast  brightness  variability.  Surprisingly,  under  some  circumstances  overcast 
twilight  2?„  exceeds  clear-sky  Ev  at  the  same  Sun  elevation.  ©  2005  Optical  Society  of  America 
OCIS  codes:  010.1290,  290.1090,  290.4210,  290.7050. 


1.  Introduction 

To  casual  observers  at  the  Earth’s  surface,  the  over¬ 
cast  sky  may  seem  to  be  the  epitome  of  visual  monot¬ 
ony,  a  nearly  featureless  shield  of  gray  cloud.  Yet 
overcasts’  transmission  of  daylight  is  seldom  gray  or 
spectrally  neutral,1  and  most  stratocumulus  over¬ 
casts  have  myriad  luminance  details  that  quickly 
move  and  change  shape  (see  Fig.  1).  Even  seemingly 
uniform  stratus  overcasts  can  exhibit  pronounced  il¬ 
luminance  changes  from  minute  to  minute,  a  fact 
long  recognized  by  lighting  engineers.2-3  More  re¬ 
cently,  researchers  analyzing  climate  and  radiative 
transfer  models  have  emphasized  that  these  models 
must  include  small-scale  spatial  and  temporal  fluc¬ 
tuations  in  overcast  optical  depth  and  transmitted 
irradiance.4-8 

Yet  despite  this  recent  work  on  overcasts’  radio¬ 
metric  variability,  we  have  found  no  comparable  re¬ 
search  that  examines  their  photometric  variability  as 
measured  from  below.  In  particular,  we  are  inter¬ 
ested  in  analyzing  short-term  fluctuations  in  overcast 
horizontal  illuminance  and  zenith  luminance  as  mea¬ 
sured  every  30  seconds  for  total  elapsed  times  tel 
<  2  h.  Our  choice  of  sampling  frequency  and  maxi¬ 
mum  measurement  period  reflects  constraints  im- 
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posed  by  (1)  our  spectroradiometers’  long  integration 
times  at  twilight  illuminance  levels  and  (2)  measure¬ 
ment  interruptions  due  to  the  onset  of  rain  or  clearing 
in  some  overcasts.  Although  our  chosen  sampling  fre¬ 
quency  would  be  too  low  to  capture  all  photometric 
details  of  the  most  rapidly  changing  partly  cloudy 
skies,9  for  many  overcasts  it  is  more  than  adequate  to 
reveal  a  remarkable  range  of  illuminance  and  lumi¬ 
nance  fluctuations. 

2.  Measuring  Overcast  Horizontal  Illuminances  and 
Zenith  Luminances 

Consistent  with  our  earlier  definition  of  “overcast,” 
we  analyze  brightness  fluctuations  only  in  those  skies 
with  complete  cloud  coverage.10  We  calculate  over¬ 
cast  horizontal  illuminances  Ev  and  zenith  lumi¬ 
nances  L„  from  our  earlier  time  series  of  >9100 
irradiance  and  radiance  spectra  measured  during  40 
different  overcasts;  integrals  are  calculated  from 
380-780  nm  in  4  nm  steps.  Our  observing  sites  are 
the  U.S.  Naval  Academy  (USNA)  at  Annapolis,  Mary¬ 
land,  the  University  of  Granada  at  Granada,  Spain, 
and  two  rural  sites  near  Owings,  Maryland  and 
Marion  Center,  Pennsylvania.  See  Ref.  1,  Sect.  2  for 
other  details  of  our  sites  and  spectroradiometers. 

Figure  2  suggests  the  wide  range  of  short-term 
variability  in  our  overcast  Ev  time  series.  In  it  we 
show  Ev  as  a  function  of  unrefracted  Sun  elevation  h0 
for  two  USNA  stratus  overcasts  on  11-21-02  and  10- 
17-03.  We  use  h0  as  Fig.  2’s  abscissa  because  for  a 
given  overcast  optical  depth  t,  overcast  Ev  at  the  sur¬ 
face  depends  on  downwelling  irradiance  atop  the 
clouds,  which  itself  is  largely  determined  by  h0.  Be¬ 
cause  we  acquired  data  after  solar  noon  on  most  days, 
we  make  h0  decrease  and  time  increase  from  left  to 
right  in  Fig.  2  and  later  figures.  For  brevity,  we  use 
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Fig.  1.  Fisheye  photographs  of  a  rapidly  changing  stratocumulus  overcast  taken  at  1750  and  1752  UTC  on  12-10-2003  at  the  U.S.  Naval 
Academy  in  Annapolis,  Maryland  (USNA).  Nikon  E5000  digital  camera  shutter  speeds  are  1/215  s  on  the  left  and  1/190  s  on  the  right; 
white  balance  and  other  exposure  parameters  (ISO  100  equivalent  sensitivity,  f/ 5.0  aperture)  are  the  same  in  both  photographs.  Light  rain 
began  falling  soon  afterward  at  —1810  UTC. 


conventional  meteorological  abbreviations  in  these 
figures:  “OVC”  for  overcast,  “CLR”  for  clear,  “St”  for 
stratus,  and  “Sc”  for  stratocumulus. 

Not  surprisingly,  the  overcast  Ev  in  Fig.  2  are  much 
smaller  than  those  for  the  9-10-03  clear  sky.  In  par¬ 
ticular,  mean  Ev  on  the  two  overcast  afternoons  is 
~3.1  times  smaller  than  mean  clear-sky  Ev  for 
27.1°  >  h0  >  21.6°.  Although  there  are  important 
exceptions  (see  Sect.  5),  the  common  observation  that 
horizontal  illuminance  is  less  beneath  overcasts  than 


Fig.  2.  Spectrally  integrated  horizontal  illuminance  Ev  as  a  func¬ 
tion  of  unrefracted  Sun  elevation  h0  for  parts  of  one  clear  and  two 
overcast  afternoons  at  USNA.  Note  that  Ev  is  scaled  logarithmi¬ 
cally. 


clear  skies  at  the  same  h0  is  a  long-lived  truism  in 
atmospheric  optics.11  Bohren  succinctly  explains  this 
commonplace  within  the  confines  of  a  two-stream 
model  of  scattering  and  absorption  by  optically  thick 
clouds.12 

Although  Fig.  2’s  two  overcasts  have  similar  mean 
Ev,  their  temporal  details  differ  greatly.  For  example, 
Ev  from  the  much  more  variable  11-21-02  overcast 
decreased  74%  in  one  19-min  period  (24.1°  >  h0 
>  22.0°),  whereas  Ev  decreased  only  —15%  on  10- 
17-03  during  the  same  h0  interval  (clear-sky  E„  de¬ 
creased  12%).  Yet  we  saw  few  visual  differences 
between  the  two  stratus  overcasts:  our  field  notes  for 
11-21-02  note  an  overcast  of  “varying  thickness  with 
some  internal  structure,”  while  that  on  10-17-03  was 
simply  “fairly  uniform.”  Neither  overcast  produced 
any  precipitation  at  the  surface. 

Beyond  merely  computing  ratios  of  consecutive  il¬ 
luminance  maxima  and  minima,  how  can  we  quantify 
Ev  temporal  variability?  One  potentially  useful  met¬ 
ric  for  comparing  the  11-21-02  and  10-17-03  overcasts 
is  D6,  the  root-mean-square  difference  between  any 
given  Ev  curve  and  a  6th-order  best-fit  polynomial  P6 
to  some  smooth  Ev  curve  (in  Fig.  3,  the  smoother 
10-17-03  F,,).13  To  determine  P6,  we  start  by  normal¬ 
izing  the  11-21-02  curve  so  that  it  has  the  same  mean 
Ev  as  the  10-17-03  curve  on  the  interval  27.0°  >  h0 
>21.5°.  Although  this  normalization  helps  reduce  D6 
for  the  slightly  darker  11-21-02  overcast,  its  D6  of 
6424  lux  is  nonetheless  36.7  times  greater  than  that 
for  10-17-03  (D6  =  175  lux).  Remarkably,  such  large 
differences  in  Ev  fluctuations  may  go  unnoticed  by  in 
situ  observers  because  overcast  Ev  often  change  over 
periods  of  minutes  rather  than  seconds. 
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Fig.  3.  Portion  of  Fig.  2’s  overcast  time  series  with  a  6th-order 
polynomial  Pe  (dashed  curve)  fitted  to  the  smoother  10-17-03  Ev 
curve.  The  root-mean-square  difference  De  between  P6  and  the 
observed  Ev  curves  is  36.7  times  greater  on  11-21-02  than  on  10- 
17-03. 


In  Fig.  3,  h0  decreases  faster  on  10-17-03  (tei 
=  35.5  min)  than  on  11-21-02  (tel  =  54.5  min;  see 
Table  1  for  more  details).  If  Fig.  3’s  abscissa  were  tel 
rather  than  h0,  the  10-17-03  time  series’  few  undula¬ 
tions  would  be  compressed  into  a  shorter  distance, 
and  thus  would  appear  somewhat  noisier.  However, 
P6  and  D6  change  only  slightly  if  we  switch  to  equal 
54.5  min  tel  that  start  at  h0  =  27.0°  on  both  days:  in 
that  case,  D6(10-17-03)  =  196  lux  and  D6(ll-21-02) 
=  5372  lux,  or  the  11-21-02  fluctuations  are  now  27.4 
times  greater.  Thus  the  11-21-02  overcast  is  much 
more  variable  whether  our  comparison  uses  h0  or  tel. 
In  fact,  we  have  determined  that  using  h0  rather  than 
tel  to  plot  overcast  time  series  does  not  affect  any  of 
our  qualitative  conclusions  here. 

3.  Fourier  Analysis  of  Overcast  Brightness  Variability 

An  obvious  weakness  in  Fig.  3’s  approach  to  quanti¬ 
fying  overcast  variability  is  that  it  does  not  preserve 
any  temporal  information  about  Ev  fluctuations.  Fur¬ 
thermore,  we  may  want  to  compare  overcasts  with 
fluctuations  that  appear  similar  in  magnitude  but 
dissimilar  in  temporal  frequency.  Fourier  analysis 
such  as  Fig.  4’s  spectral  density  function14  (SDF)  lets 
us  quantify  overcast  variability  even  if  neither  Ev 
curve  is  sufficiently  smooth  to  determine  a  meaning¬ 
ful  P6.  Note  that  hereinafter  “spectral”  means  the 
temporal  frequency  of  illuminance  or  luminance  fluc¬ 
tuations  in  a  time  series,  not  the  electromagnetic 
frequencies  in  a  single  irradiance  or  radiance  spec¬ 
trum.  SDFs  for  Fig.  2’s  overcasts  reveal  some  subtle, 
but  telling,  differences  between  them. 

In  Fig.  4,  the  influence  of  any  period  At  on  the  Ev 
curve’s  shape  increases  with  SDF,  which  itself  is  pro¬ 
portional  to  the  variance  <r2  of  Ev  at  the  given  At  (i.e., 
SDF(Af)  cr2(Ev{At))).  Note  that  Fig.  4’s  ordinate  is 


Table  1.  Sun-Elevation  and  Elapsed-Time  Intervals  for  Overcasts 


Figure 

Date 

2?„  or  L„“ 

ho  Interval  (°)6 

fel(min)6 

2 

11-21-02 

Ev 

30.7-21.6 

120.5 

2 

9-10-03 

Ev 

31.0-21.0 

53.5 

2 

10-17-03 

Ev 

27.1-21.0 

39.0 

3 

11-21-02 

Ev 

27.1-21.6 

54.5 

3 

10-17-03 

Ev 

27.1-21.5 

35.5 

5 

2-27-03 

Lv 

37.8-21.0 

120.5 

5 

3-1-03 

Lv 

38.0-21.1 

118.0 

7 

2-1-03 

Lv 

30.9-20.0 

97.5 

7 

2-27-03 

Lv 

35.0-21.0 

95.5 

9 

1-1-03 

Lv 

25.4-20.0 

59.0 

9 

2-27-03 

Lv 

26.0-21.0 

31.0 

10 

1-16-03 

Lv 

19.0-12.0 

50.5 

10 

3-30-03 

Lv 

18.7-12.3 

33.2 

11 

9-27-02 

Ev 

18.3-29.3 

61.9 

11 

10-12-03 

E„ 

22.7-18.6 

25.3 

11 

10-17-03 

E„ 

27.1-18.0 

56.5 

11 

12-10-03 

Ev 

28.1-26.4 

52.0 

12 

9-17-98 

Ev 

4.9-(— 5.5) 

53.5 

12 

10-9-98 

Ev 

4.9-(— 6.0) 

56.5 

12 

10-12-02 

Ev 

12.0-(— 5.8) 

96.5 

13 

10-10-99 

Ev 

5.1-(— 4.8) 

53.0 

13 

10-12-02 

Ev 

6.0-(-5.8) 

63.0 

13 

11-22-03 

E„ 

6.0-(-5.0) 

62.0 

14 

4-21-03 

E„ 

30.1-45.0 

80.5 

“Denotes  whether  the  given  figure  plots  illuminance  Ev  or  lumi¬ 
nance  L„. 

6For  the  given  figure,  ho  interval  and  elapsed  time  <ei  indicate  the 
period  over  which  Ev  or  Lv  is  plotted.  Other  figures  may  show  a 
slightly  different  tei  or  h0  interval  from  the  same  day.  Morning 
observations  list  the  smaller  h0  value  first,  and  numbers  in  paren¬ 
theses  denote  h0  <  0°. 


SDF  in  decibels  or  10  log10(SDF).  In  calculating 
SDF(At),  we  first  subtract  the  mean  Ev  or  Lv  from 
each  time  series,  so  the  zero-frequency  spectral  den¬ 
sity  (or  mean  power)  is  not  shown  in  Fig.  4  and  later 
SDFs. 

We  start  by  considering  in  Fig.  4  the  SDF  for  Fig. 
2’s  least  complicated  time  series,  the  monotonically 
decreasing  clear-sky  Ev  of  9-10-03.  This  SDF  is  ap¬ 
proximately  linear,  which  is  consistent  with  Fig.  2’s 
nearly  linear  decrease  in  clear-sky  Ev(h0 )  for  27.1° 
>  h0  >  2 1.5°. 15  We  say  it  is  consistent  because  a 
linear  function  such  as  Fig.  2’s  clear-sky  Ev(h0 )  can  be 
approximated  as  the  sum  of  sine  terms  whose  coeffi¬ 
cients  steadily  decrease  with  increasing  frequency  (or 
decreasing  period).  Thus  we  see  a  nearly  linear  power 
spectrum  in  Fig.  4  for  the  9-10-03  clear-sky  Ev.  Con¬ 
versely,  departures  from  linearity  in  a  SDF  (e.g.,  local 
maxima  and  minima)  indicate  that  Ev  or  L„  fluctua¬ 
tions  exist  at  the  particular  At. 

How  then  do  we  interpret  Fig.  4’s  overcast  SDFs? 
Because  the  overcast  Ev  in  Fig.  2  are  less  than  the 
clear-sky  Ev,  the  overcast  SDFs’  magnitudes  are  usu¬ 
ally  smaller.  This  happens  because  reducing  a  time 
series’  mean  Ev  tends  to  reduce  the  magnitudes  of  its 
variances  and  thus  its  SDF(Af).  However,  the  shapes 
of  Fig.  4’s  SDF(At)  curves  (as  opposed  to  their  mag¬ 
nitudes)  are  unaffected  by  a  time  series’  mean  Ev.  As 
noted  above,  local  maxima  and  minima  in  Fig.  4’s 
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Fig.  4.  Power  spectra  of  clear-sky  and  stratus  overcast  Ev  at 
USNA  (see  Fig.  2)  as  functions  of  period  A t  or  temporal  frequency 
1/At  on  the  interval  27.1°  >  h0  >  21.5°.  The  spectral  density 
function  SDF  is  defined  by  Eq.  206c  in  Ref.  14,  and 
SDF(At)  «  cr2(l?„(At))  for  variance  u2.  SDF  values  are  given  in 
dB  (10  log10(SDF)),  and  A t  is  scaled  logarithmically  to  show  more 
detail  at  smaller  periods. 


overcast  SDFs  are  caused  by  Ev  fluctuations  in  the 
corresponding  time  series.  First,  note  that  because 
the  10-17-03  SDF  has  less  spectral  detail  than  the 
11-21-02  SDF  for  At  <  10  min,  the  10-17-03  Ev  vary 
less  at  high  frequencies.  Second,  the  dB  range  of  a 
SDF  curve  at  low  frequencies  (here,  At  >  10  min) 
indicates  the  dynamic  range  or  ratio  of  maximum  to 
minimum  Ev  in  the  original  time  series.  Because  the 
10-17-03  overcast’s  Ev  dynamic  range  is  less  than 
that  for  11-21-02  (see  Fig.  3),  its  low-frequency  SDF 
range  is  also  smaller.  Third,  note  the  sparse  spectral 
detail  in  the  11-21-02  SDF  for  At  <  5  min.  This 
means  that  the  11-21-02  E„  fluctuations  with  At 
=  5-30  min  do  not  themselves  display  much  higher- 
frequency  noise. 

Not  surprisingly,  overcast  luminances  Lv  tend  to  be 
noisier  than  Ev  because  Lv  is  averaged  across  a  much 
smaller  field  of  view  (FOV):  1°  diameter  or  —2.4 
X  10  4  sr  in  our  radiometers,  compared  with  their 
2tt  sr  FOV  for  Ev.ie  Figure  5  shows  some  of  this  extra 
noise  associated  with  L„’s  narrower  FOV.  We  may 
also  be  tempted  to  infer  from  Fig.  5  that  thinner, 
brighter  overcasts  are  noisier  than  thicker,  darker 
ones:  in  the  thinner  overcasts,  the  same  change  in 
cloud  thickness  Az  should  cause  a  proportionally 
larger  change  in  cloud  transmissivity  T.12 

That  seems  to  be  true  in  Fig.  5,  where  a  brighter, 
precipitation-free  stratus  overcast  at  Owings  is  con¬ 
sistently  more  variable  than  a  darker,  snowing  one. 
But  the  appearance  of  such  time  series  can  be  deceiv¬ 
ing.  Figure  6’s  SDFs  for  these  two  skies  show  that 
their  quite  real  differences  in  temporal  L„  variability 


Fig.  5.  Spectrally  integrated  zenith  luminance  Lv  as  a  function  of 
ha  for  stratus  overcasts  at  Owings,  Maryland  on  the  afternoons  of 
2-27-03  and  3-1-03.  On  2-27-03,  light  snow  began  falling  at  the 
surface  when  h0  ~  35.7°,  and  both  the  snowfall  rate  and  snowflake 
size  increased  slightly  during  measurements.  Note  that  Ly  is  scaled 
logarithmically. 


are  less  than  we  might  think  based  solely  on  the  time 
series.  Although  the  3-1-03  overcast  is  usually 
brighter  than  the  2-27-03  overcast  (i.e.,  its  Fig.  6  SDF 
is  larger),  the  two  SDFs  are  nearly  congruent  except 
at  the  largest  and  smallest  Ah  How  do  we  interpret 
this  fact?  Since  the  SDFs  diverge  for  At  —  2-7  min 
and  the  3-1-03  SDF  has  slightly  more  spectral  detail, 
the  “extra”  Lv  variance  appears  as  high-frequency 
noise  in  the  3-1-03  time  series.  Because  the  SDFs 


Fig.  6.  Power  spectra  (SDF)  of  stratus  overcast  L„  as  functions  of 
period  A t  or  temporal  frequency  1/A t  at  Owings  on  2-27-03  and 
3-1-03  (see  Fig.  5  for  corresponding  Lv  time  series). 
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unrefracted  hQ 

Fig.  7.  Zenith  Lv  as  a  function  of  h0  for  stratus  overcasts  at  Ow- 
ings  on  the  afternoons  of  2-1-03  and  2-27-03.  Mist  and  light  drizzle 
accompanied  the  2-1-03  overcast,  and  it  has  more  luminance  noise 
than  the  snowy  2-27-03  overcast. 


converge  for  At  >  40  min,  the  time  series’  variations 
grow  more  alike  at  lower  temporal  frequencies.  For 
At  =  10-30  min,  the  nearly  congruent  SDFs  mean 
that  both  time  series  have  midfrequency  waves  of 
similar  shape. 

Figure  7  offers  a  clear  counterexample  to  our  plau¬ 
sible  idea  that  thicker  overcasts  are  less  variable. 
Here  we  compare  Fig.  5’s  snowy  stratus  of  2-27-03 
with  a  drizzling  stratus  of  2-1-03  at  Owings.  The 
latter  overcast  is  both  (1)  consistently  much  noisier 
and  (2)  alternately  darker  (thicker)  and  brighter 
(thinner)  than  the  snowy  overcast.  Perhaps  short¬ 
term  variability  depends  on  more  than  total  cloud  t: 
might  multiple  scattering  by  snowflakes  below  the 
cloud  base  give  a  snowing  overcast  smaller  Lv  varia¬ 
tions  than  a  raining  (or  precipitation-free)  overcast? 
Figure  7  suggests  that  this  idea  is  plausible. 

Figure  8’s  SDFs  for  the  Fig.  7  time  series  clearly 
show  the  high-frequency  noise  on  2-1-03  for  At 
<  5  min.  Unlike  in  Fig.  6,  at  larger  At  the  two  Fig.  8 
SDFs  have  no  matching  maxima  or  minima,  and  this 
accounts  for  the  quite  different  appearance  of  the  two 
time  series.  With  some  scrutiny,  we  can  just  divine  in 
Fig.  7’s  quite  noisy  2-1-03  time  series  the  longer- 
period  fluctuations  diagnosed  by  its  SDF  local  max¬ 
ima  at  At  =  10-12  min  and  At  ~  35-40  min.  Yet  for 
the  smoother  2-27-03  time  series,  the  lowest- 
frequency  local  maximum  in  Fig.  8  is  At  ~  30  min. 
Thus  Fig.  7’s  2-1-03  time  series  illustrates  a  more 
general  point:  high-frequency  L„  noise  may  be  super¬ 
imposed  on  a  sky  that  would  otherwise  be  less  vari¬ 
able  than  overcasts  lacking  such  noise. 

4.  Some  Physical  Explanations  of  Overcast  Variability 

Like  Fig.  7,  Fig.  9  also  shows  that  a  darker,  optically 
thicker  overcast  (1-1-03  stratus  with  fog)  is  not  nec¬ 
essarily  less  variable  than  a  brighter,  optically  thin¬ 


frequency  (min' ') 


Fig.  8.  Power  spectra  (SDF)  of  stratus  overcast  Lv  as  functions  of 
period  A t  or  temporal  frequency  1/A t  at  Owings  on  2-1-03  and 
2-27-03  (see  Fig.  7  for  corresponding  L„  time  series). 


ner  one  (2-27-03  snowing  stratus).  In  fact,  noisiness 
differences  persist  even  for  the  largest  Lv  difference 
on  the  two  days  (see  Fig.  9  at  h0  ~  22°).  Puzzling  too 
is  the  fact  that  a  fog  capable  of  obscuring  the  cloud 
base,17  surely  a  multiple  scattering  medium,  is  nois¬ 
ier  than  another  presumably  multiple-scattering  sub¬ 
cloud  medium,  the  light  snow  falling  from  the  2-27-03 
overcast. 

Yet  Fig.  10  shows  that  snow  may  have  no  consis¬ 
tent  effect  on  Lv  variability.  Here  a  nearly  visually 
featureless  stratus  overcast  on  1-16-03  is  consistently 


Fig.  9.  Zenith  Lv  as  a  function  of  h0  for  stratus  overcasts  at  Ow¬ 
ings  on  the  afternoons  of  1-1-03  and  2-27-03.  Although  the  darker 
1-1-03  overcast  looked  nearly  featureless  and  had  fog  below  its 
base,  it  also  has  more  luminance  noise. 
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Fig.  10.  Zenith  Lv  as  a  function  of  h0  for  stratus  overcasts  at 
Owings  on  the  afternoons  of  1-16-03  and  3-30-03.  Note  that  lumi¬ 
nance  noise  on  3-30-03  changes  little  when  precipitation  changes 
from  light  snow  to  light  rain. 


unrefracted 

Fig.  11.  Horizontal  illuminance  Ev  as  a  function  of  ha  for  several 
stratus  and  stratocumulus  overcasts  at  USNA  (9-27-02,  10-17-03, 
and  12-10-03)  and  at  a  rural  site  near  Marion  Center,  Pennsylva¬ 
nia  (10-12-03). 


less  noisy  than  a  precipitating  one  on  3-30-03.  Fur¬ 
thermore,  mean  overcast  t  is  similar  on  both  dates, 
and  there  is  no  obvious  change  in  the  3-30-03  time 
series’  noisiness  when  light  snow  changes  to  light 
rain  at  h0  ~  17.5°.  A  more  obvious  candidate  for 
explaining  this  noisiness  is  the  clearly  visible  cloud 
structure  on  the  latter  date. 

Does  cloud  type  (e.g.,  stratus  or  stratocumulus)  af¬ 
fect  Lv  and  Ev  variability?  Clearly  the  much  greater 
spatial  detail  of  a  stratocumulus  overcast  can  cause 
rapid  Lv  fluctuations  as  cloud  features  of  various  t 
quickly  pass  across  a  radiometer’s  narrow  FOV.  Less 
obvious  is  how  such  small-scale  features  influence  Ev 
variability,  for  which  the  radiometer  FOV  =  2tt  sr. 
Our  data  suggest  that  stratocumulus  Ev  fluctuations 
exceed  those  for  stratus  at  the  same  h0  and  cloud  t. 
But  since  these  conditions  are  rarely  met  in  nature, 
we  can  examine  only  roughly  comparable  cases. 

Figure  11  gives  four  such  cases,  and  in  them  over¬ 
cast  variability  increases  with  the  presence  of  fog 
and/or  decreasing  cloud  t.  In  our  experience, 
stratocumulus  are  often  darker  and  thicker  than 
stratus,  and  so  this  partly  offsets  the  former’s  greater 
spatial  detail.  Ultimately,  an  overcast’s  temporal 
variability  seems  to  depend  on  (and  vary  directly 
with)  its  spatial  variability.  Perhaps  surprisingly,  fog 
usually  increases  overcast  variability,  perhaps  be¬ 
cause  its  subcloud  vertical  t  fluctuates  faster  than 
does  cloud  t  itself  (assuming  that  the  fog  is  distin¬ 
guishable  from  higher-altitude  clouds).  Alternatively, 
fog  may  increase  overcast  brightness  noise  simply 
because  fog  is  closer  to  the  radiometer  than  clouds 
and  thus  fog  features’  angular  (as  opposed  to  linear) 
speed  across  the  radiometer  FOV  is  much  greater. 
This  greater  angular  speed  is  significant  because  for 
zenith  angle  0,  the  L,,(0)  contributions  to  horizontal 


Ev  °c  cos(0).  Thus  even  a  cosine-corrected  diffuser  will 
transmit  changes  in  the  angular  distribution  of  radi¬ 
ances  above  it. 

More  compellingly,  time-lapse  photographs  taken 
during  our  measurements  show  that  fluctuations  in 
overcast  L„  and  Ev  often  are  caused  by  cloud  elements 
of  different  t  that  are  advected  across  the  radiometer 
FOV  rather  than  by  t  changes  occurring  within  indi¬ 
vidual  elements.  Although  t  certainly  changes  in 
these  elements  as  they  move,  nonetheless  we  can 
often  associate  local  maxima  and  minima  in  Lv  or  Ev 
with  the  passage  overhead  of  distinct  bright  or  dark 
regions  in  the  overcast.  The  period  of  these  extrema  is 
At  <  10—15  min,  so  we  classify  them  as  mid-  to  high- 
frequency  fluctuations.  However,  the  presence  of 
small,  low-frequency  Ev  fluctuations  in  seemingly  fea¬ 
tureless  stratus  overcasts  (e.g.,  10-17-03  case  in  Figs. 
3,  4,  and  11)  indicates  that  slow  changes  in  mean 
overcast  t  also  contribute  to  overcast  variability. 

5.  Illuminance  Variability  during  Overcast  Twilights 

Still  another  plausible  idea  is  that  overcast  variabil¬ 
ity  will  decrease  during  evening  civil  twilight  because 
cloud-top  illumination  changes  from  diffuse  skylight 
plus  direct  sunlight  to  skylight  alone.  Twilight’s 
nearly  diffuse  lighting  will  reduce  Ev  differences  be¬ 
tween  lighted  and  shaded  areas  on  the  overcast’s  top, 
which  in  turn  might  reduce  Ev  fluctuations  beneath 
the  overcast.  Of  course,  as  the  Sun  nears  the  horizon, 
we  might  expect  the  ever-larger  shadows  cast  by 
nearby  clouds  to  briefly  increase  Ev  variability.  What¬ 
ever  the  merits  of  this  logic,  Fig.  12  shows  that  it 
cannot  explain  all  overcast  twilights.  Although  over¬ 
cast  Ev  fluctuations  sometimes  decrease  after  sunset 
(the  10-12-02  case),  they  also  can  either  persist 
throughout  much  of  civil  twilight  (the  10-9-98  case)  or 
scarcely  exist  before  sunset  to  begin  with  (the  9-17-98 
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Fig.  12.  Horizontal  illuminance  Ev  as  a  function  of  h0  before  and 
during  three  overcast  twilights  at  Marion  Center  (10-12-02)  and 
USNA  (9-17-98  and  10-9-98).  Nimbostratus  (Ns)  comprised  the 
base  of  the  10-9-98  overcast. 


Fig.  14.  Comparison  of  horizontal  illuminance  Ev(h0)  simulated 
by  radiative  transfer  model  MODTRAN4  and  measured  at  USNA 
on  4-21-03.  The  USNA  Ev  time  series  is  unusually  noisy,  appar¬ 
ently  because  of  persistent  thick  fog  (horizontal  visibility  distance 
—  0.5-1  km)  between  the  surface  and  the  intermittently  visible 
cloud  base. 


case).  Once  again,  angular  or  spatial  variability  in 
cloud  t  better  seems  to  explain  our  measured  Ev  fluc¬ 
tuations. 

In  fact,  Fig.  13  gives  two  spectacular  examples  of 
how  twilight  Ev  variability  depends  on  cloud  t.  In  the 
10-10-99  Marion  Center  case,  moderate  to  dense  fog 
continually  enveloped  the  observing  site,  and  pro¬ 
nounced  Ev  fluctuations  persisted  throughout  twi¬ 
light.  Note  that  fog  is  associated  once  again  with  an 


Fig.  13.  Horizontal  illuminance  Ev  as  a  function  of  h0  before  and 
during  three  overcast  twilights  at  Granada,  Spain  (11-22-03)  and 
at  Marion  Center  (10-12-02  and  10-10-99).  Near  sunset  at 
Granada,  towering  cumulus  diffusely  reflected  enough  direct  and 
diffuse  sunlight  through  the  overcast  to  briefly  make  its  Ev  exceed 
those  of  a  typical  clear  sky  at  the  same  h0  (dashed  curve). 


unusually  variable  Ev  time  series  (see  also  Figs.  9, 11, 
and  14). 

Another  extreme  is  Fig.  13’s  Granada  overcast  of 
11-22-03.  Near  sunset,  this  overcast  thinned  enough 
to  let  towering  cumulus  (Cu)  diffusely  reflect  large 
amounts  of  direct,  reddened  sunlight  to  the  surface. 
As  a  result,  the  still-intact  overcast  was  briefly 
brighter  than  a  clear  sky  at  the  same  h0  (dashed 
curve;  te\  ~  18.5  min).  This  is  the  important  exception 
noted  above  to  the  truism  that  overcasts  are  always 
darker  than  clear  skies.  In  our  experience,  such 
brightening  is  both  rare  and  visually  remarkable:  it 
occurs  only  when  tall  clouds  top  the  overcast,  and  it  is 
one  of  the  few  cases  when  overcast  Ev  fluctuations 
prompt  glances  skyward. 

At  present,  we  cannot  model  this  complex  range  of 
overcast  variability  very  realistically.  For  example, 
Fig.  14  compares  a  highly  variable  stratus  overcast  at 
USNA  on  4-21-03  with  the  Ev  trend  predicted  by  the 
radiative  transfer  model  MODTRAN4.  Despite 
MODTRAN’s  complexity  and  versatility,  its  verisi¬ 
militude  is  not  striking  in  Fig.  14  for  fixed  model 
cloud  thickness  A z.  Figure  14’s  MODTRAN  simula¬ 
tion  parameters  include  typical  seasonal  atmospheric 
profiles  of  temperature,  pressure,  humidity,  and  gas 
mixing  ratios  at  midlatitudes;  a  default  stratus  cloud 
model  from  z  =  330-930  m;  tropospheric  aerosols 
typical  of  rural  sites;  background  stratospheric  dust 
and  other  aerosols;  Mie  aerosol  phase  functions;  mul¬ 
tiple  scattering;  and  an  underlying  Lambertian  sur¬ 
face  with  reflectance  =  0.2.  Changing  any  one  of  these 
parameters  does  not  reproduce  Fig.  14’s  observed  Ev 
trend. 

Although  we  could  match  this  trend  by  continually 
varying  MODTRAN  A z,  the  exercise  would  be  purely 
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ad  hoc  without  having  detailed  knowledge  of  actual 
cloud  A z.  In  fact,  MODTRAN  makes  a  fairly  common 
modeling  assumption:  the  atmosphere  is  horizontally 
homogeneous.  Changing  a  MODTRAN  atmosphere’s 
fixed  A 2  over  time  or  h0  is  easy,  but  adding  the 
radiative-transfer  consequences  of  overcasts’  hori¬ 
zontal  inhomogeneity  will  require  new  modeling 
tracks,4-8  perhaps  including  the  use  of  fractal 
clouds.18 

6.  Conclusions 

Although  temporal  variability  of  overcast  Ev  tends  to 
increase  with  decreasing  cloud  t,  counterexamples 
are  easy  to  find  (e.g.,  Figs.  7,  9).  Twilight  h0,  subcloud 
snowfall,  and  stratus  or  stratocumulus  cloud  type 
also  have  mixed  effects  on  Ev  noisiness,  if  indeed  such 
cause-and-effect  relationships  exist.  Brightness  vari¬ 
ability  increases  more  consistently  for  (1)  fog  or  near¬ 
surface  stratus,  (2)  clouds  that  have  more  visual 
detail,  and  (3)  narrower  instrument  FOV  (i.e.,  L„  ver¬ 
sus  EJ.  Time-lapse  photography  of  overcasts  reveals 
that  their  mid-  to  high-frequency  L„  and  Ev  fluctua¬ 
tions  (say,  At  <  10-15  min)  often  are  dominated  by 
cloud  elements  of  different  t  that  move  across  the 
radiometer  FOV,  not  by  t  changes  within  individual 
clouds.  However,  lower-frequency  changes  in  mean 
overcast  t  also  appear  to  contribute  to  the  observed 
overcast  variability  (e.g.,  Figs.  3,  4,  and  11). 

The  practical  significance  of  short-term  fluctua¬ 
tions  in  overcast  luminance  and  illuminance  is  as  yet 
unknown.  However,  being  aware  of  their  existence 
and  knowing  how  to  describe  their  full  range  clearly 
is  necessary  in  developing  realistic  overcast  models. 
The  limitations  of  even  so  sophisticated  a  model  as 
MODTRAN4  (Fig.  14)  illustrate  the  large  gap  be¬ 
tween  our  ability  to  observe  and  to  predict  overcast  Lv 
or  Ev.  As  we  have  shown  above,  temporal  frequency 
analyses  such  as  the  SDF  are  useful  tools  for  quan¬ 
tifying  overcast  variability  and  for  vetting  the  time- 
and  space-dependent  details  of  new  models  of  over¬ 
cast  brightness. 
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